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Abstract 

We consider a self-gravitating system consisting of perfect fluid 
with spherical symmetry. Using the general expression of entropy 
density, we extremize the total entropy 5* under the constraint that 
the total number of particles is fixed. We show that extrema of S 
coincides precisely with the relativistic Tolman-Oppenheimer-Volkoff 
equation of hydrostatic equilibrium. Furthermore, we apply the max- 
imum entropy principle to a charged perfect fluid and derive the gen- 
eralized Tolman-Oppenheimer-Volkoff equation. Our work provides 
a strong evidence for the fundamental relationship between general 
relativity and ordinary thermodynamics. 

PACS numbers: 04.20. Cv, 04.20.Fy, 04.40. Nr 



1 Introduction 

In the past few decades, research in general relativity has suggested a very 
deep connection between gravitation and thermodynamics. The four laws 
of black hole mechanics were originally derived from the Einstein equation 
at the purely classical level[2]. The discovery of the Hawking radiation 
[5 allows a consistent interpretation of the laws of black hole mechanics 
as the ordinary laws of thermodynamics. By turning the logic around, 
Jacobson 1^ showed that the Einstein equation may be derived from the 
first law of local Rindler horizons. Inspired by Jacobson's work, a lot of 
efforts have been made to derive the dynamical equations from black hole 
thermodynamics [5]-[9]. In fact, this idea can be traced back even before 
the establishment of black hole mechanics. In 1965, Cocke [TU] proposed 
a maximum entropy principle for self-gravitating fluid spheres. Let S be 
the total entropy of spherically symmetric perfect fluid. Cocke showed 
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that the requirement that S be an extremum yields the equation of hydro- 
static equiUbrium which was originahy derived from the Einstein equation. 
However, a critical assumption in Cocke's derivation is that the fluid is in 
adiabatic motion so that the total entropy is invariant. By imposing the 
adiabatic condition, the entropy density s is expressed as the function of 
the energy density p, while for a general fluid, s is a function of at least two 
thermodynamic variables. In my opinion, the variation of S is performed 
on a spacelike hypersurface and the dynamic revolution of the fluid is ir- 
relevant. Furthermore, variation of S is not consistent with the adiabatic 
condition. If the entropy is required to be invariant, as indicated by the 
adiabatic condition, the variation of entropy would be meaningless. Thus, 
it is not appropriate and consistent to impose the adiabatic condition. In 
relation to Cocke's work, Sorkin, Wald and Zhang (SWZ)[TT] develop a 
different entropy principle for radiation. The major difference is that the 
adiabatic condition was not needed in SWZ's derivation. Moreover only 
the Einstein constraint equation was used in the proof while Cocke used 
both the constraint equation and the radial-radial component of Einstein's 
equation. However, SWZ's discussion was restricted to radiation for which 
the thermodynamic relations can be expressed explicitly and the entropy 
density only depends on one thermodynamic variable. It is important to 
know whether SWZ's treatment can be generalized to an arbitrary perfect 
fluid. In this paper, we prove a maximum entropy principle for a general 
self-gravitating perfect fluid. The new arguments used in our proof are 
as follows. First, we use the Gibbs-Duhem relation as the expression of 
entropy density for a general fluid. Second, our maximum entropy princi- 
ple is under the constraint that the total number of particles is invariant. 
Consequently, the method of Lagrange multipliers plays an important role 
in our derivation. Third, in addition to the Einstein constraint equation, 
we only make use of the ordinary thermodynamic relations to derive the 
Tolman-Oppenheimer-Volkoff (TOV) equation. No other assumptions are 
needed. Finally, we extend our treatment to a general charged fluid. With 
modified arguments, we derive the generalized TOV equation for a charged 
fluid. 



2 Review of SWZ's derivation on self-gravitating 
radiation 

Since our work is closely related to SWZ's prescription, we shall give a brief 
review on the derivation in [llj . Consider a spherical box of radiation hav- 
ing total energy M and confined within a radius R. For thermal radiation, 
the pressure p and energy density p satisfy the equation of state 

p = Ip^ (1) 
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and then the stress energy tensor is given by 

Tab = pUaUb + -pigab + UaUb) , 



(2) 



where is the 4- velocity of the local rest frame of the radiation. In terms 
of the locally measured temperature T, the energy density p and entropy 
density s are given by 



P 

s 



where 6 is a constant. So s can also be expressed as 



s = ap' 



,3/4 



(3) 
(4) 

(5) 



with a = ^b^/^. As shown by SWZ, the extrema of the total entropy S 
corresponds to a static spacetime metric 



ds^ = gtt{r)dt^ + 



2m{r) 



1 -1 



dr" + r'^dfl'^ . 



(6) 



The constraint equation, which is obtained from the time-time component 
of the Einstein equation, yields 



m'(r) 
47rr^ 



(7) 



Thus, m(r) is a mass function. 

Let the gas be confined in the region r < R. Then the total entropy is 
given by 



4-77 
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2m(r) 



2m(r) 



-1/2 



r dr 



-1/2 



r dr 



3/4 r 



2m(r) 



-1/2 



r dr . 



(8) 



Our task is to find a function m(r) such that the total entropy S is extrem- 
ized. Note that the total mass M within R is 



M = An 



r pk 

Jo 



r)r^dr = m{R) , 



(9) 



and obviously 



m(0) = . 



(10) 
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Hence, all the variations must satisfy 

Sm{0) = 5m{R) =0. 



(11) 



By using this condition, the extrema of S is equivalent to the Euler- 
Lagrange equation 



dL 
dr V dm' 



dL_ 

dm 







for the Lagrangian 

L = {m'f/^ 



2m{r) 



-1-1/2 



.1/2 



By straightforward calculation, Eq. yields 

3i. nSi. 3# 1/2 ^/ n 

m r H — m mr H — m r m r m m = U . 

16 8 8 4 2 

By substituting Eq. ([7]), one can show that Eq. is equivalent to 

d (p + p/3)[m(r)+47rr3(p/3)] 

X = i o r M ■ 

dr r[r — 2m,[r)\ 



(12) 



(13) 



(14) 



(15) 



Since p — p/3 for radiation, we see immediately that Eq. (jisp is just the 
relativistic Tolman-Oppenheimer-Volkoff equation. 



3 Maximum entropy principle for perfect fluid 

To generalize SWZ's prescription to an arbitrary perfect fluid, we first need 
to find a formula for entropy density s. Because radiation has a vanishing 
chemical potential, its entropy density depends on only one thermodynamic 
variable, e.g. T or p. For fluids consisting of particles, there are at least 
two independent variables. We start with the familiar first law 

dS = ^dE + ^dV - ^dN , (16) 

where S, E, N represent the total entropy, energy and particle number 
within the volume V. Write Eq. (fT6| in terms of densities 



d{sV) = ^d{pV) + ^dV - ^d{nV) . (17) 
By expansion, we have 

sdV + Vds = ^pdV + Vdp + j,dV - j,ndV - j;Vdn . (18) 
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Applying Eq. to a unit volume, we find 



(19) 



Combining Eqs. (fT8| and ([T9)) . we arrive at the integrated form of the 
Gibbs-Duhem relation il2, 



(20) 



To derive this formula, we only used the first law of the ordinary thermo- 
dynamics. So it is a general expression for perfect fluid. We treat (p, n) as 
two independent variables, e.g.. 



s ^ s{p,n), fi^ p.{p,n), p^p{p,n). 



(21) 



For example, the thermodynamic quantities for a monatomic ideal gas 
are given by [T3] 



—nkT , 
2 

nkT , 



P 
P 

s = —nklnT — nklnn -\ — nk 
2 2 



5 , / 2TTmk\ 

5 + '" — 



Our task is to extremize the total entropy 



s{r) 



1 



2m(r) 



-1-1/2 



(22) 
(23) 
(24) 



(25) 



In addition to the constraint Eq. (|lip . it is natural to require the total 
number of particles 



N = 4tt I n{r) 
to be invariant, i.e.. 



1 - 



2m(r) 



-1/2 



r^dr 



(57V = 0. 



(26) 



(27) 



Following the standard method of Lagrange multipliers, the equation 
of variation becomes 



SS + X5N = . 
Define the "total Lagrangian" by 



(28) 



L{m,m' ,n) = s(p{m'),n) 



1 



2m(r) 



-1/2 



Xn{r) 



2m{r) 



-1/2 



. (29) 
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Now the constrained Euler-Lagrange equation is given by 

— - 
dn ' 

d dL dL ^ 



dr dm' dm 
Thus, Eq. dSO]) yields 

ds 

— + \ = 0. 
on 

Using Eq. ([TO)) , we have 

which shows that ^ must be a constant for self-gravitating 
From Eq. (1^ . we have 

dL ( 2m\"^/^, , 



and 



Here 



dL ds o / 2m ^ 



-rM 1 - 
9m' 9to' 



ds ds dp 11 



9m' (9p dm' T Anr^ ' 
where Eqs. ([7]) and have been used. Hence 

dL 1 2m^"^/^ 



9m' 47rr V r 
and 

d (9L T{m'r — m) — r{r — 2m)T' 

dr dm' ATrT^{r — 2m)-^/^r^ 

Using Eqs. §^ and Eq. dM]) becomes 

9L / 2mY^'^fp + p 

= r 1 — 



9m \ r J \ T 
So the Eular-Lagrange Eq. ([5T|) yields 

(47rpr^ + m)r + (r - 2m)rT' = 



The constraint Eq. ([55)) yields 

// = AT' . (41) 

Rewrite Eq. (gOl) as 

p^Ts + ^in- p. (42) 

The differential of p is 

dp ~ Tds + sdT + pdn + ndp — dp . (43) 
By substituting Eq. ([T51) . we have 

dp = sdT + nd^ . (44) 
It follows immediately that 

p{r) — sT'{r) + np!(r) . (45) 
Substituting Eqs. (gSl), O and into Eq. (gS]), we have 

T' = -^p'{r) . (46) 



Substituting Eq. (j46|) into Eq. (|40|) . we obtain the desired TOV equation 

/ ^ (p + p)(47rr'^p + m) 
^ r(r-2m) ' ^ ^ 

4 Maximum entropy principle for charged fluid 

For a charged fluid, the local thermodynamic relations remain unchanged. 
For example, we can still use Eq. (l20l) as the expression of entropy den- 
sity. But the presence of charge will change the distribution of the fluid in 
spacetime. 

In coordinates (i, r, 0, 0), assume that a spherically symmetric charged 
fluid has the line element 



ds^ = gttir)df + 



^ 2m{r) ^ Q^{r) 



dr^ + r^de^ + sin^ ed(j}^ . (48) 



Here Q{r) is defined as the total charge up to the radius r and m{r) will 
be determined later. The matter field consists of a charged fluid. Let 



\/-gtt \dt 

be the four-velocity of the fluid. Then the total stress-energy tensor can 
be written as 

Tab = fab + rf,^ , (50) 
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where 



fab = pu'^u'' + p{gab + UaUb) , 



(51) 



'^ab^' ~ ^ (^a'^Pbc - -^gabF"'^Fcd 



(52) 



The electromagnetic field Fab satisfies the Maxwell's equations 

S/bF"'' = Anf = iTTpeu" , (53) 



and 



^[aFbc] =0, 



(54) 



where pe is the charge density measured by the comoving observers. By 

d 

~ dx'' 



using the identity = In-y/^ (see [Tl]) Eq. (|55| becomes 



d 



(55) 



Because of spherical symmetry, the only nonvanishing components of i^"^ 
are F*''(r) = -F''*{r). Thus, Eq. ([35]) yields 



9r(r- V-gttgrrF '') = 47rj r ^/-gttgrr = 47rpe?' 



(56) 



where Eqs. (|49)) and (l53l) have been used in the last step. So by definition, 
the function Q{r) in Eq. (|48| can be written as 



g(r) = / 4:TTr''^^/g^pedr' . 
Jo 



By comparing Eqs. ([57l) and ((56|). one finds immediately that 

1 



--Qir) 



r^V^gugr 

is a solution of Eq. (f56|) . 

Then the time-time component of Einstein's equation gives 



(57) 



(58) 



QQ' 



(59) 



This formula is consistent with the result in |15) . Now we derive the hy- 
droelectrostatic equation from the maximum entropy principle. The total 
entropy of matter takes the form 



s(r) 



2m Q2 



-1/2 



r dr . 



(60) 



8 



For simplicity, we assume all the particles have the same charge q. Thus, 
the charge density is proportional to the particle number density n 



Pe = qn. 
Together with Eq. ([57)1 . we have 



Q' 



1 



2m 
r 



21 1/2 



(61) 



(62) 



Now we treat Q{r),Q' (r) as independent variables in the Lagrangian for- 
malism. So the Lagrangian is written as 



L{m,m',Q,Q') 



1 



2m 

r 



-1/2 



(63) 



The conservation of particle number N is equivalent to the conservation of 
charge with the radius R. Now the constraints are 

m(0) = Q(0) = 0, m{R) = constant^ Q{R) ~ constant. (64) 

With these constraints, the extrema of S leads to the following Euler- 
Lagrange equations 

d dL dL 



drdQ' dQ 
d dL dL 
dm 



= 



= 



dr dm' 

To calculate the Euler-Lagrange equations, we first note that 
s = s{p, n) = s{p{m', Q, Q'), n{Q, m, Q')) . 
With the help of Eqs. ^ and (|62]), we have 

ds 
dQ' 



ds dp ds dn 



dpdQ' dndQ' 
1 Q /ill 



T Attt^ T q A'Kr^ 



2m 
r 



1/2 



Thus, 



dL 
dQ' 



1 Q 

T Aur 



^ 2m ^ Q2 



-1/2 



p. 1 1 



(65) 
(66) 

(67) 

(68) 
(69) 

(70) 



To calculate J^, first note that 

ds ds dp ds dn 

dQ ^ d^dQ^d^dQ 

1 Q' p QQ' 



T A'Kr^ T A'Kqr^ 



2m 
r 



Q' 



-1/2 



(71) 
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Then 

dL_ ^ Anr^qQsT + {fgr + ^Qfi)Q' 
dQ ~ inr^qTp/^ 

where 



(72) 



/ = !- — + (73) 



r r 

By substituting Eqs. and (O, Eq. dHS]) becomes 



= qQ^T' + Q[-mqT + qrT - qrTm + inqr^sT'^ + y^rTfiQ' -2mqrT' 
'T'] + ^r^r - 2m){fiT' - T/i') + Q^[qTQ' + ^r{fiT' - Tfi')] . 

(74) 



2r 

qr 



Using Eq. (pS)) to ehminate /i' in Eq. (|74l) . we have 

3 V7r2(r - 2m)(grr + n^^r - Tp') Q^^r[sTT' + n^^T^ - Tp') 

= qQ T -\ 1 

n n 

+ qTQ^Q' + Q[-mqT + qrT - qrTm + Anqr^sT^ + \/frTnQ' - 2mqrT' + qr^T'] . 

(75) 



Ehminating s, /i and n via Eqs. (|20l) and (|62]) . we rewrite Eq. ((75|) as 

= 47rr3(r2-2TOr + g2)(p + p)r'-47rr3(r2-2mr + g2)Tp' + rQ2Q'2 
+ QQ'{rT + 47rr3(p + p)r + g^T' + r^T' - mT - 2rmT' - rm'T) (76) 

Now we begin to calculate Eq. (l66l) . Note that 

ds ds dp 1 
dm' dp dm' Anr^T 

Then 



(77) 



dL 1 



2m g21-i/' 



dm' A-Kr^T 
Eq. dn?]) yields 

ds ds dn fi 



1 ^ ' ^2 



dm dn dm T ^^^3 _ 2m ^ 



. (78) 



(79) 
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Here we have used Eqs. and Eq. (1^ . From Eq. we find 



dL_ 

dm 



ds 
dm 



_2m Q21-1/2 



2m Q 



r + sr 
21-1 



2m Q2 

1 ^ 



-3/2 



2to g2 



2to g2 

1 ~l~ 77" 



-3/2 



-3/2 



-3/2 



r + sr 
T 47rr2(7 



iin 
T 



2to Q2 



-3/2 



1 



2to g2 



1/2 



7' 

Thus, Eq. dnH) becomes 

g^T - inr^Tip + p) + m'Tr^ - TrQQ' - rT'Q^ - r^T' 
~mrT + 2mr^T' = . 



(80) 



(81) 



Combimng Eq. (1751) and Eq. (|5T|). one can ehminate T'. Then by substi- 
tuting Eq. (15^ for to', we finally find 



gg' 

47rr^ 



{p + p) Anrp + 



TO g^^ 



2to g^ 



(82) 



This is exactly the generalized Oppenheimer-Volkoff equation for charged 
fluid [H]. 



5 Conclusions 

By applying the maximum entropy principle to a general self-gravitating 
fluid, we have derived the TOV equation of hydrostatic equilibrium. We 
only used the Einstein's constraint equation and ordinary thermodynamic 
relations. By similar assumptions but more complicated arguments, we 
have shown that the generalized TOV equation for a charged fluid can 
also be derived by extremizing the total entropy. The TOV equation is an 
important equation for self-gravitating system which was originally derived 
from the Einstein equation. Our results show that the Einstein equation 
can be derived from ordinary thermodynamic laws. This is direct evidence 
for the fundamental relationship between gravitation and thermodynamics. 



Acknowledgements 

This research was supported by NSFC Grants No. 10605006, 10975016 and 
by "The Fundamental Research Funds for the Central Universities." 



11 



References 

[I] J.D. Bekenstein, Phys. Rev. D, 7, 2333 (1973). 

[2] J.M. Bardeen, B. Carter, and S.W. Hawking, Commun. Math. Phys. 
31, 161 (1973). 

[3] S.W. Hawking, Commun.Math. Phys. 43, 199(1975). 

[4] T.Jacobson, Phys. Rev. Lett. 75, 1260 (1995) 

[5] Y.Gong, A.Wang, Phys.Rev.Lett. 99, 211301 (2007) 

[6] R.G. Cai and S.P. Kim, J.High Energy Phys. 02(2005)050. 

[7] C. Ehng, R. Guedens, and T. Jacobson, Phys. Rev. Lett. 96, 121301 
(2006). 

[8] M.Akbar and R.G. Cai, Phys. Lett. B 635, 7(2006). 

[9] R.G. Cai, and L.M.Cao, Phys. Rev. D 75, 064008 (2007) 

[10] W.J. Cocke, Ann. Inst. Henri Poincare, 2, 283 (1965). 

[II] R.D.Sorkin, R.M.Wald and Z.J.Zhang, Gen. Rel. Grav. 13, 1127 
(1981). 

[12] S.Gao, R.M. Wald, Phys.Rev. D 64, 084020 (2001) 

[13] L.D. Landau and E.M.Lifshitz, Statistical Physics (Pcrgamon, New 
York, 1981), Pt.l. 

[14] R.M. Wald, General Relativity, The University of Chicago Press, 
(Chicago, 1984). 

[15] J.D Bekenstein, Phys. Rev. D, 4, 2185 (1971). 



12 



